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SIMULATING  WATER  FLOW  IN  SOIL 
WITH  AN  ELECTRICAL  RESISTANCE  NETWORK 

By  Bunyut  S.  Vimoke  and  George  S.  Taylor2 


u 


Water  infiltration,  redistribution,  and  removal  from  soil  are  important 
processes  that  affect  all  land  management  systems.  Although  these  processes 
must  be  studied  under  field  and  laboratory  conditions,  certain  phases  can  be 
evaluated  by  analog  techniques  with  a  considerable  saving  in  time  and  cost  and 
in  accuracy.  Electrical  analogs  have  been  very  useful  in  these  evaluations. 
Two  types  have  been  used  with  good  success:  (1)  those  utilizing  the  flow  of 
electricity  through  sheets  of  electrical  conducting  paper  and  (2)  those  using  a 
network  of  electrical  resistors.  Electrical  conducting  paper  can  be  used  to 
solve  steady -state  flow  problems  in  porous  media  that  are  saturated,  isotropic, 
and  homogeneous  with  respect  to  its  hydraulic  conductivity.  An  electrical  re- 
sistance network  is  more  flexible  in  its  use  than  conducting  paper.  The  network 
can  quickly  be  adjusted  to  simulate  flow  conditions  in  soils  that  are  homogene- 
ous or  stratified,  saturated  or  unsaturated,  isotropic  or  anisotropic.  Various 
investigators  have  already  utilized  the  network  to  study  flow  of  water  either  in 
homogeneous  or  stratified  medium  (1,  6,  7,  8). 

In  this  report,  the  major  objective  is  to  present  a  method  for  calculating 
and  assembling  network  resistances  to  represent  water  flow  in  soil.  The  ap- 
proach used  herein  is  called  the  "building  block"  method.  The  soil  profile  is 
considered  to  be  composed  of  discrete,  rectangular  blocks  of  soil.  These 
blocks  are  then  joined  to  "build"  the  entire  soil  profile.  Each  block  of  soil  is 
represented  by  a  mesh  of  four  resistors.  The  meshes  are  then  joined  to  form 
the  network,  which,  in  turn,  represents  the  soil  profile.  The  discussion  in  this 
report  will  deal  with  representations  of  steady -state  flow  of  fluid  in  saturated 
and  unsaturated  media.  Procedures  will  be  given  for  representing  media  that 
are  either  homogeneous  or  stratified  and  either  isotropic  or  anisotropic  with 
respect  to  soil  hydraulic  conductivity.  Formulas  will  be  given  for  calculating 
resistances  when  rectangular  blocks  of  different  dimensions  are  used.  Special 
formulas  are  presented  for  calculating  resistances  adjacent  to  a  tubular  drain 
embedded  in  the  soil. 


Cooperative  investigations  conducted  under  cooperative  agreement  between  the  Ohio  Agricultural  Experi- 
ment Station  and  the  Soil  and  Water  Conservation  Research  Division,  Agricultural  Research  Service,  U.S.  Depart- 
ment of  Agriculture. 

Bunyut  S.  Vimoke,  formerly  Research  Assistant,  Ohio  State  University,  Columbus;  George  S.  Taylor,  Pro- 
fessor of  Agronomy,  The  Ohio  State  University,  Columbus,  and  collaborator,  Soil  and  Water  Conservation  Research 
Division,  Agricultural  Research  Service,  U.S.  Department  of  Agriculture. 

Underscored  figures  in  parentheses  refer  to  Literature  Cited  at  end  of  report. 
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CORRESPONDENCE   BETWEEN  ELECTRICAL  AND  WATER  FLOW 

IN  A  CONDUCTING  MEDIUM 

The  equations  for  the  potential  function  0  (9.)  of  fluid  flow  through  a  satu- 
rated, isotropic  medium  under  steady-state  conditions  satisfy  Laplace's  equa- 
tion, namely,  V  ^0  =  0.  Various  investigators  (4_;  <?)  have  utilized  numerical 
solutions  to  solve  Laplace's  equation  for  two  dimensional  flow  into  a  buried 
drain.  If  square  meshes  are  drawn  as  shown  in  figure  1  and  values  of  hydraulic 
head  0  are  assigned  to  each  corner,  it  can  be  shown  by  the  relaxation  method 
(9)  that  for  a  uniform  medium: 

*>0  =    <*l  +  <*2  +  ^3+^4)/4  (1) 

The  relaxation  method  can  be  applied  to  an  analogous  problem  in  electri- 
cal circuit  theory:  If  one  has  a  network  of  equal  resistances,  the  potential  VQ 
at  a  node  is  equal  to  one -fourth  the  sum  of  those  at  the  four  adjacent  nodes. 
That  is, 

v0  =  <Y,  +  y2*y3+y4)/4  (2) 

This  can  be  proved  by  utilizing  Kirchhoff's  and  Ohm's  Laws.  If  Kirchhoff's 
Law  is  applied  to  figure  2,  the  algebraic  sum  of  the  currents  i  entering  and 
leaving  the  node  of  potential  V0  is  equal  to  zero.     That  is   — 

1|     +  J2  +  i3   +  J4      =    0  ^ 


If  Ohm's  law  is  used,  equation  3  becomes   — 

(Vl  -Vp)  (V2  -Vq)        (V3-V0)     +    (V4-V0)     B    Q  (4) 

Ri  B2  B3  B4 

Solving  equation  4  for  V0>  one  obtains  the  following: 

y,/Ri  +  y2/R2  +  V3/R3  +  V4/B4  (5) 

"°  l/R,  +   l/R2+  I/B3  +    '/B4 

In  the  special  case  where  all  of  the  values  of  resistance  are  equal,  i.e.,  Ri , 
=  R2  =  R3  =  R4  =  R0  (square  network  representing  uniform  media),  the  potential 
Vo~will"T>e  given  by  equation  2.  Thus  the  voltage  V0  is  analogous  to  the  poten- 
tial 0o  used  in  the  relaxation  method. 
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FIGURE  1.- Hydraulic  head  potential 
at  corners  of  a  square  grid. 


FIGURE  2. -Portion  of  an  electrical 
resistance  network. 


GENERAL  CONCEPT  OF  A  RESISTANCE  NETWORK 

Representing  Soil  Medium  by  Electrical  Conducting  Paper 


If  one  has  a  square  piece  of  conducting  paper  with  dimension  a  by  a,  it  can 
be  used  to  represent  any  square  dimension  of  uniform  soil  medium  having  unit 
thickness  (see  appendix,  part  I).  To  illustrate,  first  select  a  piece  of  Tele- 
deltos  paper4  of  dimension  1  inch  by  1  inch.  This  piece  of  paper  can  represent 
a  block  of  soil  medium  one  foot  in  the  horizontal  direction  and  one  foot  in  the 
vertical  direction.  The  same  piece  of  paper  can  also  be  used  to  represent  a 
10  by  10-foot  block  of  soil.  In  both  cases,  the  block  of  soil  is  of  unit  thickness. 
For  example,  the  block  of  soil  represented  in  the  first  case  has  the  dimensions 
of  a  cube  of  one  foot,  and  the  block  represented  in  the  second  case  has  the  di- 
mensions of  10  by  10  feet  by  1  foot.  A  piece  of  conducting  paper  of  dimension 
5  by  5  inches  will  serve  the  same  purpose. 


Representing  Soil  Medium  by  Network  of  Resistors 


Consider  the  boundaries  formed  by  the  four  lines  that  bound  the  square  of 
resistive  paper  shown  in  figure  3a.  If  we  are  not  concerned  with  conditions 
that  exist  inside  these  boundaries,  we  would  obtain  the  same  information  from 
using  either  resistive  paper  or  a  group  of  resistors.  It  will  now  be  shown  that 
four  resistors  can  be  used  to    represent  each  square  of  resistive  paper.     Since 


A  commercial  name  for  electrical-conducting  paper.  This  type  of  paper  is  often  called  "resistive"  paper. 
The  use  of  this  or  other  commercial  products  in  this  report  does  not  imply  approval  of  the  product  to  the  exclu- 
sion of  others  that  may  also  be  suitable. 
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a  block   of   uniform  soil  can  be    represented  by    a    square  of  resistive    paper,  it 
will  also  follow  that  the  four  resistors  will  do  likewise. 
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FIGURE  3.— Ways  to  represent  a  square  piece  of  electrical  conducting  paper  by  a  group  of  resistors. 

If  groups  of  four  resistors  are  used  to  represent  a  square  block  of  soil, 
there  are  different  ways  to  put  them  together  (see  fig.  3).  The  arrangement 
that  is  most  convenient  is  shown  in  figure  3,  b,  and  this  arrangement  of  re- 
sistors will  be  used  throughout  the  rest  of  the  discussion. 

The  justification  of  using  four  resistors  to  represent  a  square  piece  of  re- 
sistive paper  can  be  shown  as  follows  (figs.  4  and  5).  Take  a  square  of  resistive 
paper  and  paint  two  highly  conductive  lines  on  it  (fig.  4).  Silver  paint  is  usually 
used,  similar  to  that  used  in  printed  circuits.  Measure  the  electrical  resist- 
ance R0  between  these  lines.  The  resistance  RQ  is  called  the  "characteristic 
resistance"  of  the  medium,  and  its  reciprocal  yields  the  electrical  conductivity 
G.  Now  do  the  same  with  the  group  of  four  resistors,  using  wires  instead  of 
painted  lines,  which  is  electrically  the  same  thing  (fig.  5).  In  all  cases,  the 
resistance  is  RQ  ohms,  which  is  the  characteristic  resistance  of  the  resistive 
paper. 
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FIGURE  4.— Square  soil  block  or  a  square  of  conducting  paper  with  two  highly  conductive  lines  painted 

on  it  to  verify  use  of  four  resistors. 
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FIGURE  5.— Square  soil  block  or  a  square  of  conducting  paper  on  which  wires  are  used  instead  of  painted 

lines  in  figure  4,  to  verify  use  of  four  resistors. 

The  directions  in  which  the  resistances  are  measured  are  orthogonal.  For 
a  two-dimensional  flow  problem,  we  can  specify  the  resistance  of  a  part  of  the 
medium  by  the  resistances  measured  in  just  two  directions.  The  resistance  of 
the  medium  in  any  other  direction  can  be  expressed  in  terms  of  its  horizontal 
and  vertical  components.  For  convenience,  one  of  these  directions  is  chosen 
to  be  the  vertical,  and  the  other ,  horizontal.  This  choice  is  made  because  many 
boundaries  in  water -flow  problems  are  usually  either  horizontal  or  vertical 
and  because  soils  are  principally  anisotropic  with  respect  to  hydraulic  con- 
ductivity in  these  two  directions. 

In  regard  to  representing  a  cross -sectional  area  of  soil  with  four  resis- 
tors, the  same  conclusion  can  be  drawn  here  as  with  the  square  of  conducting 
paper.  The  mesh  of  resistors  shown  in  figure  3,  b  can  represent  a  square 
block  of  soil  of  any  dimension.  It  is  true,  however,  that  more  information  can 
be  obtained  from  a  resistance  network  if  the  blocks  of  soil  represented  in  this 
manner  are  small.  This  is  because  one  cannot  make  electrical  measurements 
inside  the  simulated  block  of  soil. 


THE   "BUILDING  BLOCK"  APPROACH  OF  REPRESENTING  SOIL  MEDIUM 


Two  resistors  in  parallel  can  always  be  combined  into  a  single  resistor. 
The  use  of  a  single  resistor  is  merely  an  economical  convenience.  Two  re- 
sistors (Rj  and  R2)  in  parallel  are  equivalent  to  a  single  resistor  having  a 
resistance  of  R  ohms  (fig.  6). 
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B1B2 
Ri  +  B2 
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If   both   Ri    and   R2  are  equal   to    2RQ  ohms,   then   R  is  equal   to    RQ  as  shown  by 
equation  6a. 

=      B.Bz      .    2Bo-2R0     _ 
-      B,  +  B2        2Bo-2B0    "  Bo 


In  most  cases,  the  magnitude  of  a  single  resistor  in  the  network  is  derived 
from  a  consideration  of  two  resistance  values  originally  being  in  parallel.  In 
assembling  a  network  it  is  desirable  to  think  in  terms  of  the  original  resist- 
ance values  that  compose  the  building  block5  rather  than  the  values  of  single 
resistors  that  result  from  parallel  combinations.  To  illustrate  further,  look 
at  the  square  grid  ABCD  (fig.  7,  c  with  resistors  R0).  Actually,  this  is  just 
the  end  result  when  the  building  blocks  of  figure  3,  b  have  been  joined.  In  the 
assembled  network,  all  parallel  resistors  have  been  combined  into  single  re- 
sistors. It  is  important  to  remember  that  we  started  with  our  basic  building 
blbcks  as  shown  in  figure  7,  a.  If  there  had  been  no  other  resistors  in  an  ad- 
jacent mesh  to  form  a  parallel  combination,  we  would  have  left  it  unaltered; 
and  its  value  would  have  remained  at  2R0.  This  is  what  happens  at  the  bound- 
aries of  a  resistance  network  (see  fig.  8).  The  resistors  along  the  top  surface 
and  down  the  right  side  had  no  others  with  which  to  combine,  so  they  were  left 
at  their  original  value  of  2RQ.  The  reasoning  is  clear  from  the  building  block 
approach  why  boundary  resistors  are  always  twice  the  value  of  interior  resis- 
tors in  a  resistance  network. 


REPRESENTING  A  RECTANGULAR  BLOCK  OF  SOIL 


Up  to  this  time  we  have  only  considered  square  building  blocks,  these 
being  a  specialized  case  of  the  more  general  rectangular  shape.  To  determine 
the  value  of  the  four  resistors  that  represent  a  rectangular  block,  we  will  in- 
troduce a  slightly  modified  approach:  Consider  a  rectangular  block  of  soil  and 
a  piece  of  resistive  paper  of  the  same  shape.  This  block  of  soil  can  also  be 
represented  by  four  resistors  as  shown  in  figure  9.  The  problem  now  is  to 
determine  the  values  of  Ra  and  RD.  To  determine  Ra,  visualize  the  procedure 
outlined  in  figure  10.  A  piece  of  resistive  paper  is  cut  horizontally  into  two 
equal  sections.  Each  of  the  sections  of  figure  10,  b  can  be  thought  of  as  being 
represented  by  a  resistor  of  value  Ra  (figs.  10,  c,  10,  d,  and  10,  e).  This  value 
can  be  computed  by  equation  7: 

length  a  2q 

Ba  =      width      B0=     b^fi0  ■      b    *0  (?) 

The  total  resistance  in  the  horizontal  direction  is  given  by  equation  8. 

Bh  ■  Ba//Ba  ■    -g^rgj    -=t-    fBo  (8) 

5In  assembling  the  network,  the  resistors  representing  individual  blocks  of  soil  are  joined  to  form  the  entire 
network.  We  choose  to  call  this  individual  block  of  soil  and  the  mesh  of  four  resistors  that  represent  it  a  "build- 
ing block." 
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FIGURE  6.— Joining  two  resistors  in  parallel. 


(b) 

FIGURE  7. --Joining  square-grid  building  blocks. 
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FIGURE  8.— Network  of  resistors,  showing  the  value  of  the  resistance  along  the  boundaries. 
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To  obtain  the  values  for  the  vertical  resistors,  the  same  approach  is  used. 
By  taking  a  piece  of  resistive  paper  representing  our  block  of  soil  as  shown 
in  figure  11,  this  time  we  cut  it  in  the  vertical  direction.  Each  of  the  two  re- 
sulting equal  parts  can  be  thought  of  as  being  represented  by  a  resistor  of  RD 
ohms.  If  we  then  utilize  equations  9  and  10  the  total  resistance  in  the  verticaF 
direction  Ry  can  be  obtained. 

_     length  .  J±_         _    2b  (9) 


width     -o       5/2 

ffb  ffb. 
By  =     Bb^Bb    =o.  ^.D         ~      o      =nR_ 


Bb  Bb_  Bb_         t.  (io) 

Rb+Rb  2  fl^ 


To  check  the  validity  of  Ra  and  RD,  if  we  set  a  =  b,  we  then  have  Ra  -  2RQ 
and  RD  =  2R0.  The  result  agrees  with- the  resistance  of  the  square  me~sh  dis- 
cussed" previously .  The  results  of  these  calculations  are  summarized  in  figure 
12. 

We  can  now  calculate  the  proper  resistors  to  represent  a  rectangular  sec- 
tion of  soil  of  any  size.  It  is  important  to  notice  that  only  the  ratio  of  a  to  b 
occurs  in  the  final  relationships.  For  example,  the  same  resistances  ■would  be 
used  to  represent  a  block  of  soil  1  foot  by  4  feet,  2  by  8  feet,  or  3  by  12  inches. 
Each  block  has 'the  same  ratio  of  1:4.  In  figure  13  the  resistances  are  given 
for  a  number  of  commonly  encountered  examples. 


COMBINING  "BUILDING  BLOCKS"  OF  SOIL 


After  determining  the  values  of  resistances  to  represent  each  block  of 
soil,  now  we  wish  to  join  them  together.  If  we  have  a  block  of  soil  of  dimen- 
sions 2  by  2  feet  that  is  to  be  joined  to  a  2  by  4  foot  block,  the  procedure  is 
shown  in  figure  14.  To  join  side  AB  to  CD,  it  is  essential  that  AB  and  CD  must 
have  the  same  dimensions,  i.e.  of  2  units  each.  To  combine  AB  and  CD  in 
parallel,  the  resistance  to  be  used  at  this  boundary  is  given  by  equation  11. 

-ACBD     =    BAB//%) 

//  gfioBo  2   o  (ID 

=   2Bo//Bo   -    2R0+Ro  =   T fi° 

The  physical  significance  of  the  combined  blocks  shown  in  figure  14,  c  can 
be  summarized  as  follows: 

1.  The  two  blocks  now  have  an  outside  dimension  of  2  feet  by  (2  +  4  feet), 
or  2  by  6  feet. 

2.  Since  it  now  has  a  dimension  of  2  by  6  feet,  the  total  resistance  in  the 
horizontal  direction  (from  KL  to  MN)  Rh  =  3R0  and  the  total  resistance 
in  the  vertical  direction  (from  KM  to  LFT)  Ry  =  Rx>. 
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;Rb       FIGURE  9. -Representing  a  rectangular  block  of  soil  by 

four  resistors. 


(0) 


(b) 


b 
* 


-Vv*A- 


Rq 

(C) 


FIGURE    10.— Representation    of    a 

rectangular    block    of  soil   In   the 

horizontal  direction. 
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FIGURE  11.— Representation  of  a  rectangular  block  of  soil  in  the  vertical  direction. 
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FIGURE  12.— Representation  of  a  rectangular  block  of  soil  with  dimensions  a  by  b,  using  four  resistors. 
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FIGURE  13.— Representation  of  rectangular  soil  blocks  having  different  di- 
mensions with  meshes  of  resistors. 
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3.  The  outside  resistances  are  unchanged.  Each  side  portion  can  be  joined 
to  another  block  along  the  side  of  the  same  dimension.  For  example, 
KAC  can  be  joined  to  another  block  along  the  side  dimension  of  2  feet 
and  side  ACM  with  the  side  dimension  of  4  feet,  etc. 

There  is  another  way  to  determine  the  combined  resistance  AC  BD.  This 
is  done  by  a  reverse  analysis  of  the  circuit.  If  we  redraw  figure  14,  c  by  desig- 
nating resistance  AC  BD  as  Ri  (fig.  15,  a)  and  short  out  points  K,  AC,  M,  and 
likewise    points    L,  BD,  N,    we    will   have    three    resistors   joined  in   parallel   as 

shown  in  figure   15,  b.     The    overall  resistance  Ry  will  be  equal  to  =      as  given 

by  the  first  expression  of  equation  12.    We  can  obtain  the  value  of  Rl  by  utilizing 
the  remaining  expressions  in  equation  12. 

length  2  Rp 

^y  =      width     ^0   =    6    50    =    ~3~ 

p 

-y    =    2Ro^R,/R0  j12) 

or    —    +  —  +  —     =     — 
2Ro        R|        Bo  Bo 

Bi     =  -§-Bo 

To  check  the  value  of  Rh,  we  short  out  points  KL  and  MN.  We  will  then 
have  the  circuit  as  shown  in  figure  16,  a.  Notice  that  this  circuit  forms  a 
Wheatstone  bridge  so  that  the  potential  at  points  AC  and  BD  are  equal  and  we 
can  disregard  the  resistor  AC  BD.  We  now  have  the  circuit  as  shown  in  fig- 
ure  16,  b  and  which  is  analyzed  by  the  expression  in  equation  13. 


RH    =    (2R0+  4R0)//(2R0  +4R0) 


I  I  I 


"    Rh         6R0  6R0 

Bh  =  3Rq 


(13) 


The  reverse  analysis  of  a  circuit  is  very  helpful  when  we  have  to  join 
blocks  of  soil  along  sides  that  are  different  in  dimension.  The  situation  arises, 
for  example,  when  we  wish  to  expand  the  fine  meshes  of  the  resistance  network 
near  the  drain  region  to  the  coarse  meshes  in  the  region  farther  away  from  the 
drain.  This  expansion  is  done  for  both  practical  and  economical  reasons,  as 
it  has  been  found  that  in  a  drainage  system  most  of  the  potential  loss  occurs 
near  the  drain.  Since  more  detailed  information  is  needed  near  the  drain,  a 
finer  mesh  is  used  in  this  region  than  in  the  region  farther  away  from  the  drain. 

There  is  a  certain  type  of  building  block  that  will  permit  the  expansion  of 
a  fine  mesh  to  a  coarse  one.  The  building  blocks  for  this  special  boundary  are 
not  so  easy  to  simulate  as  one  may  first  think.  As  shown  in  figures  17,  a  and 
17,  b,    the    network  is    to  be    expanded  from  small   blocks    of  soil  (D  and  E)  with 
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FIGURE  14.— Joining  two  meshes  of  resistances  that  represent  two  rectangular  blocks  of  soil. 
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FIGURE  15. —"Reverse  analysis"  to  determine  the  resistance  AC  J3D.    (Also  see  fig.  14.) 
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FIGURE  16. —Reverse  analysis  to  confirm  the  horizontal  resistance  Rh  °f  tne  circuit  shown  in  figure  15,  a. 
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dimensions  1  by  1  to  larger  blocks  (F)  with  dimensions  2  by  6.  In  such  a  case, 
we  introduce  blocks,  A,  B,  and  C  into  the  picture  to  function  as  boundary  build- 
ing blocks  between  a  fine  mesh  D  or  E  and  a  coarse  mesh  F.  In  actuality, 
blocks  A,  B,  and  C  must  be  combined  to  a  single  block  with  the  right  side  to 
be  joined  to  the  coarse  mesh  and  the  left  side  to  be  joined  to  the  fine  mesh. 
The  problem  arises  as  how  to  join  a  side  dimension  of  block  A  and  B  to  the 
side  dimension  of  block  C. 


1 

D     1 

1 
A     1 

2 

C 

2 

1 
E     1 

1 
B    1 

(0) 


6R0 

■VMMr 


■VWMW- 
6R0 


(b) 


FIGURE  17.— Expansion  from  region  of  fine  meshes  to  one  of  coarse  meshes. 

To  solve  such  a  problem,  we  have  to  attack  it  indirectly.  We  cannot  say 
that  we  "join"  them  together  in  the  sense  that  blocks  were  joined  in  the  pre- 
ceding section.  Rather,  we  have  to  look  at  the  "finished  product"  and  analyze 
it  in  reverse.  For  example,  the  finished  product  of  the  boundary  building  block 
is  shown  in  figure  18,  a,  which  can  be  represented  by  a  network  of  resistors 
as  shown  in  figure  18,  b.  All  outside  resistors  are  unaltered  from  the  original 
value  of  2Ro« 

As  shown  in  figures  18,  c  and  d,  the  finished  block  has  an  outside  dimen- 
sion of  3  by  2  feet.  From  equations  8  and  10  the  total  resistance  in  the  hori- 
zontal direction  is  Rh  -  (3/2)R0  and  that  in  the  vertical  direction  is  Ry  - 
(2/3)R0.  With  the  oufside  resistances  in  figure  18,  b,  being  equal  to  2RQ,  they 
are  ready  to  be  joined  to  the  fine  mesh  on  the  left  with  side  dimension  of  1  and 
to  the  coarse  mesh  on  the  right  with  side  dimension  of  2  feet.  Now  we  must 
determine  Ri,  R2,  and  R3  in  figure  18,  b  so  that  the  overall  resistance  in  the 
horizontal  direction  Rh  is  equal  to  (3/2)R0  and  the  overall  resistance  in  the 
vertical  direction  Ry  is  equal  to  (2/3)R0,  which  are  the  original  values  of  the 
resistances  representing  a  2-  by  3-foot  block  of  soil. 
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FIGURE  18.  — Representing  boundary  building  blocks  between  fine  and  coarse  meshes  in  the  resist- 
ance network. 

In  determining  Ri,  R2,  and  R3,  we  have  only  two  equations,  i.e.,  the  equa- 
tions for  Rh  and  Ry.  However,  it  is  possible  to  set  Rl  =  R3  because  of  sym- 
metry. We~~now  have  the  network  as  shown  in  figure  19,  a.  To  formulate  the 
equation  for  Ry,  we  short  out  points  AEB  and  DFC.  We  will  then  have  the  net- 
work as  shown  in  figure  19,  b.  Notice  also  that  the  circuit  G,  AEB,  H,  DFC 
forms  a  Wheatstone  bridge  so  that  the  potential  at  points  G  and  H  are  equal. 
Therefore,  we  can  disregard  the  resistance  R2  in  the  calculation.  We  now 
have  the  circuit  as  shown  in  figures  19,  c  and  the  magnitude  of  R]  is  obtained 
by  the  following: 


or 


(2R0  +   2R0)//(R|  +  R|)//2R0   =    B, 


I 


4R, 


2R 


Ql 


2R 


no 


3 
2R0 


B\  -   3  Bo 


3  -0 


(14) 


To  solve  for  R2,  short  out  points  AGP  and  BC  of  figure  19,  a.  The  new 
network  is  shown  in  figure  20,  a.  This  network  can  also  be  thought  of  as  two 
sets  of    equivalent    resistors  AEB    and   AFB,  connected  in   parallel  (See  fig.  20, 
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FIGURE  19.  — Circuit  analysis  to  determine  the  Resistance  Ri  at  the  boundary  between  fine  and  coarse  meshes  in  a 

network. 
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FIGURE  20. —  Circuit  analysis  to  determine  the  resistance  R2  as  shown  in  figure  18,  b  and  in  figure  19,  a. 
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b  and  c),    so   that  each  set   will  have    a    resistance  of  2Rf}.     Therefore,  set  AEB 
can  be  analyzed  by  equation  15  (fig.  20,  d  and  20,  e): 


[(2R2+     -§-R0)//2R0]    +   2R0   ■  2Ry  ■  3R0  (15) 

Notice  that  until  now  nothing  is  said  about  the  inside  dimensions  of  GH, 
EH,  and  HF  of  figure   18,  a. 

An  interesting  point  is  observed  with  regard  to  node  H  in  the  interior  of 
this  building  block  (fig.  21,  a).  The  question  arises  as  to  where  to  plot  the 
value  of  the  voltage  measured  at  this  internal  node.  One's  first  inclination 
might  be  to  plot  it  at  the  intersection  designated  by  H  in  figure  21,  b.  It  can 
be  shown,  however,  that  the  actual  point  of  correspondence  is  located  at  X. 
Consider  the  situation  shown  by  figure  21,  c.  Let  us  short  out  AGD  and  BC  and 
then  apply  3  volts  at  BC  while  AGD  is  attached  to  the  ground  (i.e.,  at  zero  volt). 
If  we  measure  the  potential  at  E  or  F,  we  will  have  a  reading  of  1  volt,  indi- 
cating that  both  E  and  F  are  located  one  third  of  the  distance  from  AD  to  BC 
as  expected.  If  we  then  measure  the  voltage  at  H,we  will  have  a  reading  of  2/3 
volt,  indicating  that  GH  (or  in  reality  GX  of  fig.  21,  b)  is  2/3  the  distance  of 
either  AE  or  DF.  The  potential  of  2/3  volt  at  point  H  can  also  be  calculated, 
if  we  use  Kirchhoff's  law  having  the  voltages  at  adjacent  nodes  E,  F,  and  G  as 
1,  1,  and  0  volt,  respectively,  with  the  resistance  EH,  FH,  and  GH  each  being 
(2/3)RQ.  Therefore,  in  the  analog  system,  the  voltage  measured  at  the  node 
H  in  figure  21,  a  should  be  plotted  at  point  X  located  2/3  the  distance  between 
GH,  as  shown  in  figure  21,  b.  The  finished  product  of  the  2-  by  3-foot  block  of 
soil  is  shown  in  figure  21,  d,  and  its  network  of  resistors  is  shown  in  figure 
21,  a.  This  same  circuit  was  described  by  Liebmann  (6_,  pp.  97,  figs.  10  and 
and  11).  If  one  analyzes  Liebmann' s  circuit,  one  will  find  that  Rv  -  (2/3)R0. 
However,  Rj-j  of  Liebmann' s  circuit  will  be  equal  to  (5/3)Ro  instead  of  (3/2)RQ, 
which  represents  a  discrepancy  of  11.1  percent.  Also,  if  a  potential  of  3  volts 
is  applied  across  BC  and  AGD  of  Liebmann' s  circuit,  one  would  find  that  the 
voltage  at  point  E  or  F  is  1.2  volts  and  the  voltage  at  point  H_  is  0.6  volt.  These 
values  indicate  that  the  distance  AE  or  DF  is  2/5  of  the  distance  AB  instead  of 
the  value  of  1/3. 

The  value  of  the  resistance  GH  of  figure  21,  a,  can  be  used  to  locate  point 
X  of  figure  21,  d.  In  all  cases  the  resistance  RAE  1S  equal  to  RDF-  The  dis- 
tance GX  is  related  to  distance  AE  as  shown  in  equation  16. 

GX    =    (2Rgh/Rae)AE  (16) 

Further    illustrations    will   be    shown   later    on   in  the  discussion   of   anisotropic 
soil. 
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FIGURE  21.— Locating  the  interior  point  X  in  the  soil  block  that  corresponds  to  the  point  H  in  the  network. 


With  the  information  on  different  types  of  building  blocks  discussed  above, 
we  are  now  ready  to  put  them  together  to  form  a  network  representing  a  ho- 
mogeneous soil  medium.  For  example  let  us  represent  a  section  of  soil  having 
dimensions  of  5  by  20  by  1.  The  dimensional  units  are  arbitrary,  but  we  will 
consider  the  soil  to  have  dimensions  of  5  by  20  feet  by  1  foot.  The  soil  may  be 
represented  by  the  building  blocks  of  soil  shown  in  figure  22,  a  with  side  di- 
mensions as  indicated.  Figure  22,  b  gives  the  resistances  representing  the 
building  blocks,  while  figure  22,  c  shows  the  assembled  network.  The  arrange- 
ment of  resistance  shown  in  figure  22,  c  will  yield  the  most  accurate  informa- 
tion on  potentials  along  the  left-hand  boundary  where  the  meshes  are  smallest. 
This  conclusion  follows,  as  potentials  in  a  network  can  only  be  measured  at  the 
nodes.  This  arrangement  is  particularly  useful  when  a  subsurface  drain  is 
located  in  the  region  of  small  meshes.  Since  the  potential  drop  is  quite  large 
near  a  sink  (i.e.,  the  drain),  greater  accuracy  is  needed  in  that  region. 

In  a  later  section  (p.  24)  the  simulation  of  building  blocks  in  the  region 
near  a  circular  drain  will  be  discussed.  It  will  be  shown  that  those  resistances 
which  represent  the  drain  can  be  inserted  at  any  depth  on  the  left  side  of  the 
network  shown  in  figure  22,  c.  If  a  half-drain  were  placed  on  the  left  boundary 
and  constant  (but  different)  potentials  applied  at  the  drain  and  at  the  top  bound- 
ary of  the  network,  one  would  simulate  ponded  flow  into  a  drain.  The  left-hand 
side  of  the  network  would  represent  the  vertical  plane  through  the  center  of  the 
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FIGURE  22.— Representing  homogeneous  soil  sections  with  a  resistance 

network. 

drain.  The  top  part  of  the  network  would  represent  the  soil  surface,  and  the 
bottom  part,  an  impervious  layer.  The  right  boundary  of  the  network  would 
represent  a  vertical  plane  across  which  flow  does  not  occur  and  would  corre- 
spond to  the  plane  between  equally  spaced  drains  in  level  topography  (9). 
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Observe  that  in  figure  22,  a,  if  we  use  a  mesh  size  of  1  foot,  we  will  have 
the  impervious  layer  at  5  feet.  If  we  wish  to  locate  the  impervious  layer  at  4 
feet,  we  can  remove  the  top  row  of  building  blocks.  With  this  type  of  arrange- 
ment, we  can  build  the  network  with  odd  or  even  depth  to  the  impervious  layer 
by  adding  or  removing  the  top  row  accordingly. 


REPRESENTING  A  STRATIFIED  SOIL 


So  far,  all  the  concepts  deal  with  homogeneous  medium  where  the  entire 
soil  is  of  the  same  hydraulic  conductivity  K.  In  actuality,  the  soil  profile  is 
divided  into  different  horizons.  Generally,  each  horizon  has  a  different  con- 
ductivity. For  simplicity,  we  will  choose  a  soil  profile  having  three  layers  of 
different  conductivities:  Kl,  K2,  and  K3,  respectively,  from  the  top  layer  to 
the  bottom  one.  Also,  we  will  designate  Kl,  K2,  and  K3  as  2.24,  0.56,  and  0.28 
cm./hr.  The  ratios  Kl  :  K2  :  K3  will  be  8:2:1.  These  particular  values  are 
chosen  so  that  K1/K2  =  4  and  K2/K3  =  2.     Let  the  characteristic  resistance  R0 


K,  =    2.24 


(a) 


nWW- 


FIGURE  23. —Representing  square  blocks  of  a  stratified  soil  having  hydraulic 
conductivities  Ki,  K2,  and  K3  with  meshes  of  resistors. 
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be  5,000  ohms  and  equal  to  R2  m  order  to  represent  K2.  We  can  see  that  Ri 
(representing  Ki)  is  5,000/4  or  1,250  ohms  and  R3  (representing  K3)  is  5,000 
x  2  or  10,000  ohms.  In  the  Kl  region,  a  square  block  of  soil  can  be  represented 
by  four  resistors  as  shown  in  figure  23,  b.  Figures  23,  e  and  g,  represent  a 
square  block  of  soil  in  the  K2  and  K3  regions,  respectively.  Since  5,000  ohms 
is  used  as  characteristic  resistance  of  a  square  mesh,  it  is  more  convenient 
to  express  both  Ri  and  R3  in  terms  of  R£.  Figures  23,  c  and  h,  show  the  trans- 
formations. An  example  of  the  networks  for  stratified  layers  is  shown  in  fig- 
ure 24. 

A  more  general  expression  for  calculating  Rh  and  Ry  is  given  in  the  ap- 
pendix, part  II,  for  the  case  where  rectangular  meshes  of" different  sizes  are 
employed  and  the  hydraulic  conductivity  also  differs  in  the  various  horizontal 
layers. 


REPRESENTING  AN  ANISOTROPIC  SOIL 


In   the   previous    discussion  we    were  dealing    with  soil  that   had   a  different 
hydraulic    conductivity    in   each   layer    or    horizon.     It   is    sometimes   found  that 
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(d) 


FIGURE  25.— Representing  square  blocks  of  soil  having  different  hydraulic  conduc- 
tivities in  the  vertical  and  horizontal  directions  (anisotropic  soil). 
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FIGURE  26.— Representing  sections  of  ani- 
sotropic soil  having  two  layers  with  a 
resistance  network. 
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within  the  same  layer  of  soil,  the  hydraulic  conductivities  in  the  horizontal  and 
in  the  vertical  directions  are  also  different.  This  type  of  soil  is  called  aniso- 
tropic with  respect  to  its  hydraulic  conductivity.  To  set  up  a  network  to  rep- 
resent an  anisotropic  soil,  let  us  designate  a  soil  with  two  layers.  The  upper 
layer  has  a  hydraulic  conductivity  in  the  vertical  direction,  Kyi  as  0.12  cm./hr., 
and  in  the  horizontal  direction,  KhI  as  0.03  cm./hr.  In  thelower  layer,  Kh2  = 
0.06  cm./hr.  and  Ky2  =  0.03  cm."7hr.  By  designating  the  conductivity  of  "0.03 
cm./hr.  as  1,  we  can  see  that  the  ratio  Ky:KH  of  the  upper  layer  is  4:1  and 
of  the  lower  layer  is   1:2  (fig.  25).    By  choosing  Fhe  characteristic  resistance  of 
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5,000  ohms  (R0)  to  represent  the  conductivity  of  KH1  and  KV2  (each  equaling 
0.03  cm./hr.),we  can  see  that  resistance  Rvi  representing  Kyi  is  1,250  ohms, 
whereas  the  resistance  RH2  representing  ~KH2  is  2,500  ohms.  In  the  upper 
layer,  a  square  block  of  soil  can  be  represenTed  by  four  resistors  as  shown  in 
figures  25,  a  and  b,and  in  the  lower  layer  by  figures  25,  c  and  d.  An  assembled 
network  of  resistors  for  this  soil  is  shown  in  figures  26,  a,  b,  and  c.  .The  over- 
all distances  of  the  simulated  flow  medium  are  arbitrarily  selected. 

Special  attention  is  called  to  the  location  of  the  internal  nodes  of  the  bound- 
ary building  block  between  the  coarse  mesh.  As  shown  in  figure  26,  a,  the  in- 
ternal node  of  the  boundary  building  block  of  the  upper  layer  is  plotted  at  the 
point  11/12  of  1  and  that  of  the  lower  layer  at  1/3  of  1.  These  locations  must 
be  determined  by  either  calculation  or  experimentation  as  shown  in  a  previous 
section,  p.   16. 
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FIGURE  27. —Diagram  showing  the  approach  used  to  calculate  the  resistance 
R2  see  (e)  which  represents  the  soil  region  around  a  drain.  The  value  £2  is 
equal  to  8R,  where  R  is  the  resistance  of  the  conducting  paper  shown  in  (b). 
The  value  of  R  is  evaluated  by  using  the  analogy  between  R  and  the  imped- 
ance Z0  of  the  transmission  line  shown  in(c).  The  relationship  between  R_ 
and  Zo  is  expressed  by  equations  17  and  18. 
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REPRESENTING  A  CIRCULAR  DRAIN 


Up  to  this  point  building  blocks  always  represent  rectangular  sections  of 
soil.  In  studies  dealing  with  subsurface  drains,  however,  the  drain  is  usually 
circular  in  cross  section.  To  cover  this  case,  we  must  extend  the  building 
blocks  concept  to  include  such  curved  sections.  To  approach  this  problem, 
first  cut  from  a  piece  of  resistive  paper  the  section  drawn  with  a  solid  line  as 
shown  in  figure  27,  a.  In  this  figure,  r  is  the  drain  radius  and  s  is  the  square 
mesh  size  around  the  drain.  At  the  present  time  it  will  be  convenient  to  re- 
member that  this  section  is  just  one  fourth  of  the  entire  section  as  shown  in 
figure  27,  a.  If  we  paint  the  boundaries  of  the  square  and  the  circle  with  con- 
ducting paint  as  shown  in  figure  27,  b,  we  could  then  measure  the  resistance 
R.  However,  direct  measurement  is  sometimes  neither  convenient  nor  accu- 
rate. This  problem  has  been  solved  mathematically  for  electrical  transmis- 
sion lines  and  wave  guides.  The  resistance  R  of  figure  27,  b,  is  essentially  the 
same  as  the  impedance  Z0  of  a  square  wave  guide  with  a  circular  inner  con- 
ductor as  shown  in  figure  27,  c.  The  only  difference  between  R  and  Z0  is  a 
conversion  factor  Zq  relating  the  characteristic  resistance  of  free  space  and 
the  characteristic  resistance  (resistance  per  square  =  RQ)  of  the  resistive 
paper.     The  relationship  is  given  by  equation  17  (5,  pp.  426-429). 

R    -    —  Z  <17) 


The  expression  yielding  Z0  for    a   transmission  line  as  shown  in  figure  27,  c  is 
given  by  equation  18  (12,  p.  590). 

Z0*    138    \oq\oP    +  6.48  -  2.34A  -  0.48B  -  0.I2C 

where   />  =  -§-     (Fig.  27b) 

I   +  0.405/?"4 
I    -  0. 405/0 ~4 

I   +  O.I63/>-8 
I    -  O.I63/>-8 


A    = 


B    = 


(18) 


r    .      I    +  0067  />-'2 
-   "      I   -  0.067/3-12 

and   Zq    =  Jp      =**    I20tt   =  376.7  ohms    where 

P-0    =   permeability   of   free  space 
Sq     -   permitivity   of  free    space 
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In   our    network,    if   we    are    dealing   with   a    square  mesh   around  the    drain 

with  mesh  size  s_  and  radius  of  the  drain  r,we  see  that  s_  =  =  and  r_  =  =.    There - 

D       s 
fore  p  =  ^p  =  =.     From   these  expressions,    we    can  determine  R  of   figure  27,  b. 

d       jr.  r  —  ° 

Returning  to  the  original  problem,  if  we  paint  conducting  lines  on  the  out- 
side and  inside  boundaries  of  figure  27,  d  and  then  measure  the  resistance,  we 
will  get  4R  ohms  since  this  part  represents  only  one  of  the  four  quarters  that 
are  in  parallel,  as  shown  in  figure  27,  a.  The  resistance  of  the  section  shown 
in  figure  27,  d  can  be  represented  by  the  resistance  shown  in  figure  27,  e.  In 
this  particular  measurement,  the  resistors  of  value  Ri  are  effectively  shorted 
out  so  that  the  circuit  reduces  to  figure  27,  f.  Taking  an  overall  view  of  figure 
27,  d  leads  us  to  the  relationship: 

R2 
4R  *     9i  or 

ezo  <19) 

R2   =   8R   =     2{j    Bo 


The  next  question  to  ask  is  what  value  to  pick  for  the  Ri  resistors  in  fig- 
ure 27,  e?  In  our  studies,  the  value  of  Rl  =  2RQ  is  chosen.  In  every  way  this 
value  seems  reasonable,  but  there  is  no  basis  to  assume  that  it  is  exact.  It 
may  very  well  depend  on  the  diameter  of  the  drain  as  does  R2.  However,  a 
number  of  experiments  using  R\  equal  to  2RQ  have  yielded  results  (11)  which 
are  in  good  agreement  with  exact  analytical  solutions  (3,  equation  11). 

To  illustrate  the  above  approach  in  calculating  the  resistances  to  repre- 
sent the  drain  region,  let  us  designate  the  size  of  the  mesh  as  1  foot,  radius 
of  the  drain  r_  =  3  inches,  Ro  =  5,000  ohms. 

Zo     =   138  log,^    +  6.48  -2.34A  -  0.48B  -  0.I2C 

12 
where   p    =  -5-    =  4 

A  =  1.00317 
B  »  I 

C%  I  <19a> 

Zo    =    138  logio4  +  6.48  -  2.3474  -  0.48  -0.12     =    86.6169 

r       =   |4oRo  =     8X86.6169      x   500Q    =  9|9?  Qhms 
_c  Zq    —  376.7 


2 


=    4599  ohms 
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From   the    above    example,    we    can    see  that   for    a  given   value    of  p   we  can 

8ZQ  8Zq 

calculate  the  drain  resistance  R2  =  — j —  Ro-    If  we  define  the  constant  Cd  =  , 

-      'Z0     -  --      Z0 

the  drain  resistance  is  given  by  R2  =  CdRo«  Table  1  reports  a  few  values  of 
Cd  along  with  the  corresponding  values  of  p  .  The  values  reported  therein  will 
suffice  for  calculating  those  values  of  the  drain  resistance  that  are  normally- 
encountered  in  network  studies. 


Table  1 .  -Some  values  of  the  constant  Cd  as  a  function  of  p.  The  function 
p  is  equal  to  the  ratio  ^/_r,  where  £  is  the  size  of  the  square  mesh  at  the  drain 
and  r  is  the  drain  radius.  (See  fig.  27,  a.)  The  magnitude  of  the  drain  resist- 
ance R2,  as  shown  in  figure  28  is  given  by  R2  =  CdR0 


s 

P   -  8-2o 

s 

8Z0 
C  d  =  ZcT 

1.0000 

0.00319 

2 . 1818 

1.06632 

1.0435 

0.07491 

2.4000 

1.18823 

1.0909 

0.14598 

2.6667 

1.32278 

1.1429 

0.21300 

3 . 0000 

1.47298 

1.2000 

0.28227 

3 . 4286 

1.64316 

1.2632 

0.35321 

4.0000 

1.83949 

1.3333 

0.42635 

4.8000 

2.07162 

1.4118 

0.50261 

6.0000 

2.35568 

1.5000 

0.58247 

8.0000 

2.72187 

1.6000 

0.66677 

9.0000 

2.88452 

1.7143 

0.75628 

12.0000 

3.23794 

1.8462 

0.85195 

24.0000 

4.12018 

2.0000 

0.95483 

00 

00 

We  now  have  the  building  blocks  necessary  to  represent  the  area  around 
the  drain.  In  figure  28  we  put  these  building  blocks  together  to  form  a  network 
in  a  homogeneous  soil  with  a  square  mesh  around  the  drain.  (See  also  Appendix, 
pt.  II.) 
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FIGURE  28. —  Diagram  showing  connection  of  resistance  R2   which  represents  the  drain 
region  to  those  representing  square  soil  blocks. 


AN  EXAMPLE  UTILIZING  THE  RESISTANCE  NETWORK 
TO  STUDY  FLOW  PROBLEMS 


In  network  studies,  there  are  three  things  in  which  we  are  primarily  in- 
terested: (1)  the  potential  distribution  from  which  we  can  obtain  —  among 
other  things  —  equipotential  lines;  (2)  the  streamlines  from  which  we  can  learn 
the  direction  of  water  movement  through  the  soil;  and  (3)  the  flow  rate  or,  in 
the  case  of  drainage,  the  amount  of  water  that  we  can  remove  from  the  soil 
during  a  given  period  of  time.  All  these  data  can  be  obtained  directly  with  the 
network. 

As  an  example,  consider  the  ponded  flow  case  in  drainage:  Multiple  drains 
of  radius  r  are  equally  spaced  in  a  homogeneous  soil  at  a  distance  a  and  at  a 
depth  d.  Water  is  ponded  on  the  level  ground  surface  to  a  height  t.  The  drains 
are  completely  permeable  to  water  and  are  running  full  with  no  back  pressure. 
An  impermeable  horizontal  layer  lies  at  a  depth  h  below  the  ground  surface. 
First,  the  network  is  assembled  as  shown  in  figure  29.  To  obtain  the  potential 
distribution,  an  e.m.f.  (electromotive  force)  of  100  volts  is  then  applied  to  the 
top  boundary.  The  drain  terminal  (fig.  28)  is  connected  to  the  ground  as  shown 
by  the    diagram   in  figure  29.     A   voltmeter  is  used  to    measure    the  potential  at 
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each  node.  A  vacuum  tube  voltmeter  (VTVM)  is  recommended  since  its  in- 
ternal resistance  is  very  high  (around  1 1  megohms)  and  the  potential  measure- 
ment can  be  made  accurately.  For  convenience,  the  potential  difference  of  100 
volts  is  applied  to  the  upper  network  boundary  so  that  the  voltage  reading  at 
any  node  yields  the  potential  directly  as  a  percentage  of  that  at  the  ground 
surface.  After  the  potentials  are  measured  at  each  node,  the  equipotential  lines 
can  be  drawn  by  interpolating  between  potentials  at  adjacent  nodes. 

To  obtain  the  flowline,  the  boundaries  of  the  network  are  now  "reversed" 
as  shown  in  figure  30.  This  is  done  by  first  removing  the  attachments  from 
the  top  boundary  and  the  drain  terminal.  This  time  an  e.m.f.  is  applied  to  the 
previously  unconnected  boundaries  as  shown  in  figure  30.  The  potentials  are 
again  measured  at  each  node,  and  a  new  set  of  equipotential  lines  are  drawn. 
The  new  set  of  equipotential  lines  will  be  the  flow  lines;  that  is,  the  two  sets  of 
equipotential  lines  will  be  orthogonal. 

To  measure  the  drain  flow  rate,  two  different  procedures  can  be  followed: 
From  equation  44  in  appendix,  pt.  I,  we  see  that  the  applied  voltage  V  and  the 
hydraulic  head  potential  0  are  related  by  a  conversion  coefficient  C0-  The 
potential  at  the  top  boundary  0n  is  equal  to  d  +  t,  where  d^  is  the  vertical  dis- 
tance between  the  drain  center- and  the  soil  surface  and  t  is  the  thickness  of 
ponded  water  on  the  surface.  The  potential  at  the  drain  0O  is  equal  to  r,  the 
radius  of  the  drain.  This  boundary  condition  is  obtained  from  the  assumption 
that  the  drain  is  running  full  with  no  back  pressure.  Thus,  the  potential  dif- 
ference from  the  top  boundary  to  the  drain  is  equal  to  0n  -  0O  or  d  +  _t  -  r_.  If 
we  apply  V  volts  to  the  upper  network  boundary,  the  potential  Vn  at  the  top 
boundary  is  V.  The  potential  at  the  drain  terminal  Vp  is  zero.  Thus  (Vn  -  Vo) 
=  V  volts;  and  CQ  can  be  determined  by  equation  Z0. 


Co     ■ 


d  ♦  t  -r 


(20) 


FIGURE  29.— Application  of  voltage  to  the  network.  The 
positive  pole  of  the  battery  is  connected  to  the  top 
boundary  and  the  ground  Is  connected  to  the  drain.  A 
vacuum  tube  voltmeter  (Y)  Is  used  to  measure  voltage 
distribution  at  each  grid  point. 
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FIGURE  30.— Reversing  the  boundary  of  the  network  to 
obtain  streamlines. 

After  Co  is  determined,  we  then  measure  the  total  current  I  passing  through 
the  circuit.  The  flow  rate  Q  is  then  determined  by  equation  48  in  appendix,  pt. 
I. 

One  disadvantage  of  the  above  method  is  the  error  and  inconvenience  in- 
volved in  measuring  the  current  I.  In  the  first  place,  the  characteristic  re- 
sistance Ro  should  be  chosen  in  the  range  of  5,000  to  10,000  ohms  to  minimize 
errors  in  voltage  readings  at  the  network  nodes.  (See  also  appendix,  pt.  III.) 
This  range  is  chosen  so  that  the  resistance  of  the  network  will  be  far  less  than 
the  internal  resistance  of  the  VTVM  of  1 1  megohms  and  more  than  the  re- 
sistance due  to  wiring  and  metal  contact  (i.e.,  from  1  to  10  ohms).  When  the 
characteristic  resistance  is  chosen  in  this  range,  the  current  passing  through 
the  network  will  be  too  low  for  accurate  measurement.  To  increase  the  current 
reading,  one  has  to  increase  the  applied  voltage.  This  is  somewhat  dangerous 
and  inconvenient.  Also  there  is  the  risk  of  altering  the  resistance  values  in 
the  network  due  to  greater  current  flow. 

These  disadvantages  are  avoided  by  using  the  following  method  for  deter- 
mining the  flow  rate:  the  total  resistance  Rn  between  the  top  boundary  of  the 
network  and  the  drain  terminal  is  measured"  This  is  done  by  attaching  wires 
A  and  B  of  figure  29  to  an  ohmmeter  rather  than  to  a  battery,  or  better,  to  a 
Wheatstone  bridge.  The  resistance  of  the  network  Rn  can  be  measured  to  an 
accuracy  within  0.5  percent  with  a  Wheatstone  bridgeT  The  current  I_  is  given 
by   V    ;  ancj  from  equation  20,  V  =  Co  (d+t-r).      Therefore,   the    current  I  can  be 

expressed  as  follows: 


C0(d  +  t-  r) 


Rn 


(21) 
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If  we  substitute    the   value  of  I  in  equation  48   of   appendix  I  and  designate    Q1  as 
the  flow  rate  of  the  whole  drain,  we  obtain  the  relationship  given  by  equation  22. 

0'     .   2K(d  +  l  -l)  n  (22) 

9  fin  -° 

The  factor  2  shown  in  equation  22  is  to  convert  the  flow  rate  from  that  of 
a  half  drain  to  that  of  a  whole  drain  Q'.  Also  from  equation  22  we  can  see  that 
the  drain  flow  rate  Q'  is  independent  of  the  applied  voltage  V  and  the  current  I. 
All  we  need  to  know  is  the  resistance  Rn  across  the  network.  The  general 
expression  for  the  flow  rate  can  be  expressed  as 

0     .     K(0a-»O)    gQ  (") 

Rn 

where  0n  is  the  hydraulic    potential  at  the  top   boundary    and  0O  is  the  hydraulic 
potentiarat  the  drain  or  at  the  point  of  discharge. 

The  drawings  shown  in  figures  29  and  30  are  only  schematic.  For  more 
detail,  let  us  assemble  a  network  to  represent  a  flow  problem  for  the  following 
dimensions:  The  half  spacing  between  drains  (a/2)  is  72  feet,  and  an  impervi- 
ous layer  (h)  is  18  feet  below  the  ground  surface.  The  depth  to  the  drain  center 
line  (d)  is  12  feet.  The  drain  diameter  (2r)  is  4  inches.  A  water  table  is  main- 
tained at  the  ground  surface,  and  a  hydraulic  conductivity  (K)  of  0.001  foot/min. 
is  assumed. 

Because  of  symmetry,  the  network  is  assembled  to  represent  only  half  of 
the  flow  region.  For  horizontal  distances  up  to  6  feet  from  the  drain,  each 
linear  foot  in  the  soil  is  represented  by  one  resistive  unit  (i.e.,  a  1  -  by  1 -foot 
mesh).  Between  6  and  8  feet,  a  2-  by  2-foot  mesh  is  used,  while  a  2-  by  8-foot 
mesh  is  employed  for  distances  greater  than  8  feet.  For  a  2-  by  8-foot  mesh, 
the  horizontal  distance  is  8  feet.  A  potential  difference  of  100  volts  is  applied 
through  the  network,  the  potential  at  each  node  is  measured,  and  the  equipoten- 
tial  lines  are  evaluated  by  interpolating  the  potential  between  adjacent  nodes 
(see  fig.  31).  After  reversing  the  boundaries  and  again  measuring  the  poten- 
tials, the  flow  lines  are  determined  in  a  similar  manner  and  are  also  shown  in 
figure  31. 

The  characteristic  resistance  R0  is  chosen  to  be  5,000  ohms,  and  the  re- 
sistance across  the  network  Rn  for  the  drain  depth  at  12  feet  is  measured  to 
be  8,797  ohms.    Using  equation~22,  we  determine  the  flow  rate  Q  to  be 

.  2  X  0.001   (12+0  -  2/12)  X  5,000 
=  8,797 

=  0.01345  (ft.3/ft.  drain-min.) 
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In  the  case  of  stratified  or  anisotropic  medium,  the  flow  rate  Q  can  be 
determined  the  same  way  as  for  homogeneous  medium.  However,  we  have  to 
keep  in  mind  that  the  value  for  the  hydraulic  conductivity  K  is  that  chosen  to 
be  represented  by  the  characteristic  resistance  RQ.  For  example,  consider 
the  case  for  flow  of  water  through  a  stratified  soil  having  hydraulic  conductivi- 
ties of  2.24,  0.56  and  0.28  cm./hr.,  respectively  (see  section  VI).  If  R^  is 
5,000  ohms  and  represents  the  conductivity  K^  of  0.56  cm./hr.,  the  values  of 
K.2  and  R2  will  appear  in  place  of  K  and  R0,  respectively,  in  equation  22.  The 
value  of  Ki  and  K3  can  be  ignored,  since  their  overall  effect  is  already  in- 
cluded in  the  resistance  across  the  network  Rn.  In  figure  32  a  flow  net  for  the 
flow  of  water  through  an  anisotropic  soil  is  shown.  The  method  for  calculating 
the  resistances  is  given  in  the  section  beginning  on  p.  21. 

The  results  obtained  from  an  electrical  resistance  network  are  associated 
with  errors  that  are  due  to  the  mesh  size  employed,  the  selection  of  the  char- 
acteristic resistance  R0,  and  the  inaccuracies  in  setting  individual  resistors. 
The  errors  due  to  mesh  size  is  discussed  by  Liebmann  (6),  whereas  those 
errors  caused  by  the  other  two  sources  are  discussed  in  the  appendix,  part  III. 


REPRESENTING  AN  UNSATURATED  REGION  OF  SOIL 


Under  saturated  flow  conditions ,  the  hydraulic  conductivity  K  in  each  layer 
of  soil  is  independent  of  its  location  within  the  layer.  In  unsaturated  flow,  the 
hydraulic  conductivity  of  a  particular  layer  may  differ  from  place  to  place. 
For   flow   problems    encountered   during    drainage,   unsaturated   flow    generally 
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FIGURE  31.— Flow  net  for  a  homogeneous,  Isotropic  soil  during  ponded  flow  Into  a  4- inch  drain. 


31   - 


T 


5J-r 
I 


4^1 


STREAM  I 


4K 


POTENTIAL 


2K 


■*§6 


4  6  8 

Distance  from  center  of   drain   (ft.) 


10 


12 


FIGURE  32.— Flow  net  for  a  layered,  anisotropic  soil  during  ponded  flow  into  a  4-inch  drain. 

occurs    above  a   nearly   horizontal   water   table  whereas    saturated  flow   occurs 
below  this  surface. 

To  represent  a  block  of  soil  having  variable  spatial  conductivity  with  a 
mesh  of  resistors,  first  consider  a  piece  of  conducting  paper  with  the  side 
dimension  a  by  b  as  shown  in  figure  33.  Cartesian  coordinates  are  set  up  with 
the  x-axis  along  the  lower  edge  and  the  y-axis  along  the  left  boundary.  Let  us 
assume  that  the  resistivity  R  is  only  a  function  of  the  vertical  coordinate,  i.e., 
R  =  R(y).  This  condition  approximates  that  of  a  small  block  of  unsaturated  soil 
in  which  changes  in  moisture  content  (and  thus  in  hydraulic  conductivity)  in  the 
horizontal  direction  are  quite  small,  as  compared  to  those  in  the  vertical. 
Previous  work  has  shown  this  condition  to  be  essentially  the  case  for  unsatu- 
rated flow  above  a  falling  water  table. 

If  we  use  the  "building  block"  approach,  we  must  calculate  the  overall 
resistance  Ry  in  the  vertical  direction  and  the  overall  resistance  Rh  in  the 
horizontal  diFection.  Since  the  resistivity  varies  only  in  the  y_  direcTion,  let 
us  then  consider  a  small  strip  with  dimension  dy_  by  a_.  The  incremental  re- 
sistance ARy  can  be  expressed  as 


ARy  =    dy  R(y) 
a       — 


(24) 


From  inspection,    one  can  see  that  the  incremental    resistances  ARy  are  joined 
in  series,    so   that  the  total    resistance    in   the  vertical  direction  Rv~is    equal  to 
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SARy.     Therefore, 


Ry  ■  ±.J  V-y 


(25) 


To  determine  Rh>    consider  the  same    strip   in  figure  33.     The  incremental 
resistance  ARfj  is  expressed  as 


ARH  =  ■£-  B(y) 
dy      - 


(26) 


Since   the   total    resistance    Rg   is    the    result   of  joining    in   parallel   the    incre- 
mental strips  of  resistance  ARh>    direct  integration  of  equation  26  is    not  valid. 


FIGURE  33.— Diagram  showing  flow  region  with  variable  resistivity  in  the  ^-direction.  The  re- 
sistivity is  symbolized  by  R(y),  with  values  at  the  upper  and  lower  boundaries  given  by  Rj  and  Ri, 
respectively. 
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However,    if  we  designate  the  termAGH  as  the    incremental   conductance  in  the 

horizontal  direction  and  equal  to  L_  equation  26  can  then  be  expressed  as 

ARH 


so  that 


AGH=  i^L  (27) 

_      a  R 


Gu   =_L 


^H 


Bh 


(28) 


From  the  procedure  given  above,  two  points  should  be  emphasized:  (1) 
Integration  of  incremental  resistances  can  be  made  if  they  are  joined  in  series; 
(2)  integration  of  incremental  conductances  can  be  made  if  they  are  joined  in 
parallel. 

To  illustrate  the  use  of  equations  25  and  28,  first  assume  that  the  resistiv- 
ity R(v)  varies  linearly  in  the  y_  direction  between  y_  =  O  and  y_  =  b .  This  as- 
sumption introduces  only  a  small  error,  unless  the  function  R(y)  is  changing 
rapidly  in  the  vicinity  O  <  y_  <  b  or  the  soil  block  represented  is  of  relatively 
large  dimension,  i.e.,  a  coarse  mesh  is  used.     Thus, 

R  =    R(y)=  M  +  B  (29) 

where  A  and  B  are  constants.  From  figure  33,  designate  the  resistivity  of  the 
lower  edge  as  Ri  and  of  the  upper  edge  as  R j .  (The  determination  of  Ri  and  Rj 
will  be  discussed"  in  a  later  section.)  By  seTting  first^  =  o  and  ^  =  b,  the  con~ 
stants  A  and  B  can  be  determined  as 

A  =  E-i_-  Bl  (30) 

b. 

B  =   Ri  (3D 


Therefore, 


(32) 
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Substituting  R(y)  of  equation  32  in  equations  25  and  28  and  carrying  out  the 
integrations,  we  fina  that 


Ry  =  b    Ri±Ji 


(33) 


Bh=^ 


a(Bj_-Bj) 

£  loge  Si 
Ri 


(34) 


To  test  the  validity  of  equations  33  and  34,  let  us  consider  the  case  when 
the  entire  flow  region  is  saturated.  In  such  case  Ri  and  Rj  will  be  equal  to  Ro, 
the  characteristic  resistance  of  the  square  piece  of  resistive  paper  having 
uniform    resistivity  R0.     We   then  need  to  prove   that  under  such  conditions  Ry 

~       b  a  ~ — 

and  Rh  also  be  equal  to-^^R0  and-^-R0>  respectively.  It  is  obvious  from  equa- 
tion 33  that  if  Ri  and  _R_i  are  replaced  by  R0thenR"V  =  —  Rp .  The  same  substitu- 
tion, however,  will  not  yield  direct  proof  in  the  case  of  Rh*  For  if  Rj  and  Ri 
are  replaced  by  Ro  in  equation  34,  we  would  have  the  indeterminate  form  _. 
With  the  help  of  equations  30  and  31,  equation  34  can  be  rewritten  as 


Bh  = 


aA 


loge 


iB 


(35) 


If  Rj  =  Ri  =  R0,  then  from  equations  30  and  31  we  would  have  A  =  0_  and  B  =  RQ. 
If  the    values  of   A  and  B    are    substituted   into    equation  35,    Rh   will  be  equal  to 

rr-;  and  L'Hospital'  s  Rule  is  applicable  with  respect  to  A.  Therefore,  differen- 
tiations of  the  numerator  and  the  denominator  of  equation  35  with  respect  to 
A  yields 


Rh  = 


Mm 
A-*0 


aA 


loge 


(JiA  +  H 
B 


Mm 
A-^0 


B 


=  ai=£B0  (36) 


It  should   be    noticed  that   the    expression   Rh  =-1F  B    in   equation  36  is    valid 
both  in  the  saturated  and  unsaturated  zones,  as  long  as  Rj   =  Ri  =  B.  Accordingly, 
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in  the    saturated  zone    B    is  equal   to  R0,   whereas  in   the   unsaturated  zone  B  is 
equal  to  either  Rj  or  Ri. 

After  Ry  and  Rh  are  determined,  we  can  now  proceed  to  represent  the 
block  of  soiFor  resistive  paper  shown  in  figure  33  with  four  resistors.  This 
representation  is  shown  in  figures  34,  a  and  b.  The  terms  Ry  and  Rh  are  cal- 
culated from  equations  33  and  34,  respectively.  —  ~ 


(a) 


(b) 


FIGURE  34.— Representing  the  flow  region  shown  in  (a)  with  a  mesh  of  4  resistors. 

From  equations  33  and  34,  which  yield  the  values  of  Ry  and  Rh,  we  have 
two  unknowns  to  determine  —  Ri  and  Rj.  These  terms  represent  the~resistivity 
of  the  lower  and  upper  boundaries,  respectively,  of  the  nonuniform  conducting 
paper  shown  in  figure  33.  To  evaluate  Ri  and  Rj  for  the  study  of  flow  in  un- 
saturated soil,  it  is  necessary  to  relate  the  resistivity  of  the  conducting  paper 
R  and  the  hydraulic  conductivity  of  the  soil  K. 

In  unsaturated  flow,  the  hydraulic  conductivity  is  dependent  on  the  mois- 
ture content  of  the  medium.  As  the  moisture  content  decreases  below  satura- 
tion, the  hydraulic  conductivity  likewise  decreases.  Also,  as  the  moisture 
content  decreases,  the  hydrostatic  pressure  h  of  the  soil  water  films  increases 
in  the  negative  direction.  In  other  words,  the  soil  moisture  tension  is  in- 
creased. The  unsaturated  hydraulic  conductivity  is,  in  turn,  related  to  the 
hydrostatic  pressure.  The  relationship  between  conductivity  and  hydrostatic 
pressure  can  be  obtained  by  methods  described  by  Childs  and  Collis -George 
(2_)  and  by  Neilsen  and  coworkers  (10).  An  example6  of  this  type  of  relationship 
is  shown  in  figure  35.  From  the  graph  one  can  see  that  each  type  of  medium 
has  a  characteristic  relationship  between  the  hydraulic  conductivity  and  the 
hydrostatic  pressure. 

To  relate  the  resistivity  R  to  the  hydrostatic  pressure  h,  we  first  use  the 
relationship  between  the  resistivity  R  of  the  conducting  paper  and  the  hydraulic 

c 

conductivity  K  of  the  fluid  flow  medium.     This  can  be  written  as  K=-«f-  where    £ 


6Thiel,  T.  J.    Unpublished  Master's  Thesis,  Ohio  State  University,  Columbus,  Ohio  1959. 
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FIGURE  35,-Fluid  conductivities  K  of 
fine  and  coarse  sands  as  a  function  of 
hydrostatic  pressure  h.  The  points 
plotted  on  and  to  the  right  of  the  or- 
dinate represent  the  saturated  fluid 
conductivities  of  the  two  sand  sizes. 


is    a    constant  relating  K  to  R.     From  this  relationship  we  can  then  write  equa- 
tion 37. 


lo.i 
R    '  Ko 


(37) 


K0  and  K  are  the  hydraulic  conductivities,  respectively,  of  saturated  and  un- 
saturated soil;  and  RQ  and  R  are  the  electrical  resistivities,  respectively,  of 
uniform  (analogous  to  saturated  soil)  and  nonuniform  conducting  medium  (anal- 
ogous to  unsaturated  soil).  Both  K  and  R  are  space  functions  in  their  respective 
media.  Another  source  of  needed  information  is  the  equation  relating  the  po- 
tential head  0,the  elevation  head  y_,  and  the  hydrostatic  pressure  head  h.  These 
three  are  related  by  equation  38  for  cases  where  only  hydrostatic  pressure  and 
elevation  components  of  potential  are  significant. 


0  =  y  +  h 


(38) 


With  the  information  above,  we  are  now  able  to  utilize  the  network  analog 
to  solve  unsaturated  steady-state  flow  problems.  To  illustrate,  consider  a 
soil  section  having  a  dimension  of  10  by  10  feet,    as  shown  in  figure  36.    On  the 
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right  and  left  side  of  the  soil  section,  the  hydraulic  heads  are  kept  constant  at 
8  and  6  feet,  respectively.7  For  simplicity ,  assume  that  the  characteristic  re- 
lationship between  K/KQ  and  the  hydrostatic  pressure  h  is  as  shown  in  figure 
37.  From  figure  37,  we  have  K/K0  =  1  for  h  >  -  1.0,"^/Ko  =  0.5h  +  1.5  for 
-2.8  <_h  <  -1 ,  and  K/K0  =  0.1  for  h  _<  -2.8.  The  same  problem  can  be  assembled 
on  the  network  with  each  mesh  of  resistors  representing  a  small  square  block 
of  soil  of  dimension  1  by  1  foot.  The  grid  points  are  shown  in  figure  36.  Let 
us  first  assemble  the  network  to  represent  a  saturated  and  homogeneous  soil 
medium  having  resistances  R0  everywhere  in  the  interior  and  2Ro  along  the 
boundaries.  Eight  boundary  resistors  on  the  right  side  and  six  on  the  left  side 
are  shorted  out  as  shown  in  the  figure,  and  potentials  of  8  and  6  volts  are  ap- 
plied along  the  right  and  left  boundaries,  respectively.  The  potential  at  each 
grid  point  is  measured  and  recorded.  The  magnitude  of  h  at  each  grid  point 
is    then   calculated   by    equation    38.     From   the    values    of   h,    Ri     and   Rj  can  be 


'A  similar  type  of  study  had  been  reported  by  Luthin  and  Day  (8). 
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FIGURE  36.— Network  representing  a  10'  x  10'  section  of  soil.  Each  square  mesh  has  a 
dimension  of  1'  x  1'.  Hydraulic  head  potentials  of  6  and  8  feet  are  kept  along  the  left 
and  right  boundaries,  respectively.  These  boundary  potentials  are  obtained  by  applying 
an  e.m.f.  of  6  and  8  volts  along  the  left  and  right  boundaries  of  the  network. 
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FIGURE  37.— Illustrative  curve,  showing  the  relationship  of  K/Ko  and  hydrostatic  pressure  h.    The  saturated 

conductivity  is  given  by  Ko. 

obtained  from  the  relationship  given  in  figure  37.  The  resistances  of  the  un- 
saturated medium  can  then  be  calculated  as  illustrated  by  the  diagram  shown 
in  figure  38.  The  section  of  soil  represented  in  figure  38  is  that  shown  with  a 
heavy  line  in  the  interior  of  figure  36. 

To  calculate  the  resistances,  consider  the  section  in  figure  38,  a,  with  the 
potentials  0  recorded  at  each  node.  As  explained  above,  the  initial  values  of 
the  potential  are  obtained  by  utilizing  a  homogeneous  and  saturated  flow  me- 
dium. The  next  step  is  to  adjust  the  resistors  in  accordance  with  the  magni- 
tude of  the  hydrostatic  pressure  h.  To  proceed,  the  section  in  figure  38,  a  is 
first  separated  into  four  small  squares  (see  fig.  38,  b).  The  potential  at  each 
corner  is  also  shown  in  figure  38,  b.  By  knowing  the  potential  0  and  the  ele- 
vation from  the  datum  plane  y_,  the  hydrostatic  pressure  can  be  calculated. 
For  example,  the  potential  at  the  upper  left  corner  of  block  No.  1  is  6.86  feet 
and  y_  is  9  feet.  Therefore,  the  hydrostatic  pressure  h  is  6.86  -  9  feet  =  -2. 14 
feet  as  shown  in  figure  38,  c.  From  the  known  value  of  h,  the  ratio  of  con- 
ductivity K/K0  can  be  obtained  from  figure  37  and  is  found  to  be  0.43.  The 
same  procedure  is  applicable  for  obtaining  K/Ko  when  the  experimental  curve 
shown   in   figure  35  is  used   in   the  network    study.     However,    if   the  problem  is 
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solved  by  numerical  analysis  and  a  high  speed  electronic  computer  is  used,  the 
relationship  between  K/Ko  and  h  must  be  expressed  in  an  analytical  equation, 
the  latter  sometimes  being  somewhat  cumbersome  to  use.  This  difficulty  can 
be  avoided  by  approximating  the  curve  shown  in  figure  35  with  a  series  of  broken 
straight  lines  for  which  linear  equation  can  be  assigned  to  each  line. 

After  the  value  of  K/KQ  is  obtained,  the  value  of  R/R0  can  be  calculated 
from  equation  37  as  1/0.43  =  2.33  and  is  shown  in  figure  38,  d.  With  the  same 
procedure,  the  value  of  R/Ro  at  the  upper  right  corner  of  block  No.  1  is  calcu- 
lated to  be  2.04  (figure  38,  d),  and  the  value  of  the  two  are  averaged  to  yield 
2.19.  This  value  is  then  set  equal  to  Rjl/Eo  so  that  R  j  1  =  2.19  Ro,  as  shown  in 
figure  38,  e.  The  magnitude  of  Rii  of  block  No.  1  can  D~e  calculated  in  the  same 
manner  starting  from  the  values  of  0  at  the  lower  corners.  From  Ril  and  Rj  1 , 
the  values  of  Rvi  and  RH1  oi  block  No.  1  can  be  calculated  with  equations*33 
and  34  and  are  found  to  be  1.62  Ro  and  1.55  Ro,  respectively ,  as  shown  in  figure 
38,  f.  The  network  of  four  resistors  representing  block  No.  1  is  shown  in  figure 
38,  g.  The  same  procedure  is  carried  out  for  blocks  Nos.  2,  3,  and  4,  and  they 
are  joined  together  in  parallel,  as  shown  in  figure  38,  h.  The  same  procedure 
is  carried  on  throughout  the  unsaturated  region  of  soil. 
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FIGURE  38.  — Diagrams  illustrating  calculation  of  the  resistances  which  repre- 
sent unsaturated  flow  medium.    See  text  for  detailed  discussion. 
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FIGURE  39.— Flow  net  for  steady-state  flow  of  water.  Water  is  ponded  at 
left  and  right  boundaries,  and  the  relationship  between  _K  and  h  is  that  shown 
in  figure  37.  The  procedure  illustrated  in  figure  36  and  figure  38  is  followed 
on  the  analysis. 


At  this  point  the  voltages  are  then  reapplied  to  the  network,  and  the  poten- 
tial distribution  is  again  recorded.  The  new  values  of  node  potential  are  shown 
in  figure  38,  i.  One  can  see  that  the  potential  distribution  of  the  new  set  is  dif- 
ferent from  the  previous  one.  The  entire  procedure  is  repeated  until  no  sig- 
nificant potential  change  is  observed,  and  the  relaxation  is  completed.  The 
flow  net  for  this  particular  problem  is  shown  in  figure  39. 

The  discussion  above  deals  only  with  an  unsaturated  medium  whose  con- 
ductivity varies  only  in  the  vertical  direction.  If,  in  some  cases,  the  conduc- 
tivity varies  in  both  vertical  and  horizontal  directions,  the  determination  of 
Ry  and  Rfj  can  only  be  obtained  approximately.  It  would  appear,  however,  that 
acceptable"  accuracy  can  be  obtained  by  using  successive  linear  approximations 
of  resistivity  R.  On  this  basis  the  expressions  Ry  and  Rh  for  a  block  of  soil 
having  dimensions  a  by  b  can  be  approximated  as  — 


*.*§+& 


(39) 


41   - 


and 


!.h-  tm 


where  Ri  and  Rj  are  determined  similar  to  those  illustrated  in  figure  38.  This 
approximation  is  based  on  the  assumption  that  the  unsaturated  soil  is  homoge- 
neous with  respect  to  the  conductivities  in  the  horizontal  and  vertical  direction. 
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APPENDIX 


I.    REPRESENTING  WATER  FLOW  IN  SOILS  WITH  ELECTRICAL 
CURRENT  FLOW  IN  CONDUCTING   PAPER 


A  square  of  conducting  paper  with  dimensions  a  by  a  can  be  used  to  rep- 
resent any  square  dimension  of  uniform  soil  which  has  unit  thickness.  For 
example,  consider  the  block  of  saturated  and  homogeneous  soil  of  hydraulic 
conductivity  K  shown  in  figure  40,  a.  Assume  that  constant  potentials  0n  and 
0O,  respectively,    are    applied  to  the  upper    and   lower  surfaces.    A  steady  "state 
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FIGURE  40.  — Diagram  illustrating  the  analogy  of  water  flow  through  soil  and 
electrical  flow  by  conducting  paper.  One- dimensional  flow  is  illustrated  in 
(a)  and  (b);  and  2-dimensional  flow  in  (c)  and  (d). 
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flow  will  be  established,  and  the  quantity  of  water  Q  leaving  the  lower  surface 
per  unit  time  will  be  given  by  equation  41. 

Q  =  y   A    !=    "K^-     A  (41) 

where  v  is  the  average  flow  velocity  per  unit  area  and  A  is  the  dimension  of 
the  square  block  in  the  x-direction.  Similarly,  the  current  I_  leaving  the  lower 
edge  of  the  conducting  paper  shown  in  figure  40,  b  will  be  given  by  the  follow- 
ing relationship: 

I      =      U       =        "SfyT         °  (42) 

where  G  is  the  electrical  conductivity,  i  is  the  current  flux  per  unit  length  in 
the  x'-direction,  and  V  is  the  voltage. 

Since   both   media  are  homogeneous    with    respect  to  their  conductivity,  the 
functions  0  and  V  will  be  related  in  the  following  manner: 


c0(^jti-  *j)  =<yj+i-Yj) 


(43) 


where  Co  is  a  proportionality  constant.  By  evaluating  equation  43  for  (J_  +  1)  = 
n  and  j  =  0,  it  can  be  seen  that  the  voltage  V  is  related  to  0  by  equations  44  and 
45. 

C0    =    (Yn  -Yo)/(*n-*o) 

(44) 


A 

(45) 


where  _?'=  I    -jr 


i  _       -fl. 
1  - 1  X 


C,     =    (Yo  -C0^o) 

For  the  examples  shown  in  figures  40,  a  and  b,  it  can  be  seen  by  inspection 
that  both  0  and  V  vary  only  in  the  vertical  direction.  Thus,  if  equation  45  is 
differentiated  with  respect  to  y_,  the  following  results: 

5/  \hl    '  Qo  dy  (46) 
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By  substituting  equation  46  in  equation  42  and  then  dividing  the  resulting  ex- 
pression into  equation  41,  one  obtains  the  following  relationships  among  the 
conductivities  and  flow  rates  in  the  two  media: 


Q   = 


LK  (47) 

C0G 


Q  -    -=*-,  <48> 

I 

where    R0  =   ~q 

From  the  above  analysis,  it  can  be  seen  that  a  square  of  electrical  con- 
ducting paper  can  represent  one -dimensional  flow  of  water  in  any  square  di- 
mension of  uniform  soil  having  unit  thickness.  The  potential  0  is  related  to 
the  voltage  V  as  shown  in  equation  45,  and  the  flow  rate  Q  is  given  in  terms  of 
G,  K,  and  I  by  equations  47  and  48.  Although  the  example  presented  here  is 
for  one -dimensional  flow,  equation  45  and  47  are  equally  valid  for  two-dimen- 
sional problems.     Examples  of  the  latter    are  illustrated  in  figures  40,  c  and  d. 

An  approach  similar  to  that  used  to  derive  equations  45  and  47  can  be 
followed  to  show  that  any  rectangular  piece  of  conducting  paper  of  dimension 
a  by  b  can  represent  any  rectangular  block  of  soil  of  dimension  A  by  B.  Ex- 
pressions similar  to  equations  45  and  47  can  thus  be  derived  which  contain  the 
parameters  a,  b,  A  and  B. 


II.    STEPS  TO  DETERMINE  THE   POTENTIAL  AT  A  NODE  EST 
NUMERICAL  ANALYSIS  OR  IN  A  RESISTANCE  NETWORK 

One  of  the  more  diverse  cases  encountered  in  determining  the  potential 
at  a  node  in  saturated,  isotropic  soil  occurs  when  the  node  is  (1)  located  at  a 
point  on  the  interface  of  horizontal  layers  with  different  hydraulic  conductivity 
K  and  (2)  is  surrounded  by  "building  blocks"  of  different  sizes.  Figure  41 
shows  the  location  of  such  a  node.  From  figure  41,  a,  four  blocks  of  soils  have 
dimensions  a  by  b,  b  by  c_,  a  by  d,  and  d  by  c,  respectively.  The  two  blocks  at 
the  top  have  a  conductivity  Kj  and  the  two  blocks  at  the  bottom  have  a  con- 
ductivity Kfi.  The  ratio  of  Kt7K.b  is  given  by  m.  If  R0  is  the  characteristic 
resistance  ~bf  the  upper  blocks,  the  four  blocks  can  be  represented  by  four 
meshes  of  resistances  as  shown  in  figure  41,  b.  The  parallel  sides  of  the  mesh 
can  be  joined  together  as  shown  in  figure  41,  c  and  the  final  combination  is 
shown  in  figure  41,  d.  In  figure  41,  d,  Ri  is  the  combined  resistance  of  RAB 
and  RCD,  R2  is  that  of  RDM  and  RGN>  etc.  If  Vl  is  the  potential  at  node  AC, 
V2  is  the  potential  at  node  MN,  etc.,  then  VQ  is  the  potential  at  node  BDGE. 
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FIGURE  41.— Building  block  approach  used  to  calculate  the  potential  _V0  in  relation  to  the  voltages  at 
adjacent  nodes,  the  sizes  of  rectangular  blocks  surrounding  the  node  of  potential  V0,  and  hydraulic  con- 
ductivities of  the  blocks.    (See  text  in  appendix,  part  II,  for  discussion.) 

The    potential  Vo  is   now   to  be  determined.     Using  Kirchhoff '  s  law,  we  have 
the  following: 


Vi  -Vo  v2  -Vo    .      V3-V0    +     V4-V0 

Ei  B2  B3  B4 


(49) 


Now  substitute  Rl,  Rz,  R3,  R4  into  equation  49  and  solve  for  Vo.     We  then  have 


since 


2b  Ro  2mc_Ro 

-1  (a  +  c)'  -2     =     fmb  +  d) 


R3 


2mdRo 


R4 


2moRo 


Yo  = 


(a  +  c)  '       ~4  (mb  +  d)  (50> 

mocd(a  +  c)yi  +  obd(mb  +  d)\/2  +  obc(o  +  c)V3  +  bcd(mb  +  d)\/4 
macd  (q  +  c)    +    abd(mb  +  d)    +    qbc(a+c)    +   bed  (mb  +  d) 
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If  V0  is  within  the  layer  where  Kj  -  Kg    so    that  m  =  1 ,  and  also  if  a  =  c_,  b  =  d, 
and  a/b  =  r,  we  will  have  the  following~expression: 

(y,  ♦  y3)   »  (72)  (V2  ^4)  (51) 

v0  ■ 


0  +  ^) 


which  agrees  with  results  previously  reported  elsewhere  (4_,  equation  8,  p.   39). 
Finally,  if  a  =  b  so  that  r  =  1,  equation  51   reduces  to  equation  52 

V|    +  N/2    +V3  +  V4  (52) 

Yo  ■ 4 > 

which  is  the  same  as  equation  2. 

The  use  of  Kirchhoff's  law  to  determine  the  potential  at  a  node  as  de- 
scribed above  is  also  applicable  in  the  case  of  unsaturated  steady  state  flow. 
As  an  example,  let  us  first  refer  to  figure  38:  By  designating  the  potential  of 
6.95  volts  at  the  center  node  of  the  network  fig.  38,  i,  as  V0>  one  will  find  that 
this  value  can  also  be  obtained  by  Kirchhoff's  law.  This  is  done  by  designating 
the  four  surrounding  potentials  of  6.97,  7.12,  6.93,  and  6.78  volts  as  Vl,  Vz> 
V3,  and  V4,  respectively,  and  likewise  by  designating  the  resistances  of  1.53 
R0,  1.17  R0,  1.02  R0,  and  1.24  RQ  as  Rl,  R2,  R3,  and  R4,  respectively  (see  fig. 
38,  h).  By  substituting  these  values  into  equation  49  and  solving  for  V0,  one 
will  find  that  the  value  Vo  is  6.95  volts  which  is  the  same  as  the  voltage  ob- 
tained by  measurement. 

Potential  at  the  node  adjacent  to  the  drain 

In  general,  the  building  blocks  around  the  drain  region  are  square.8  Let 
us  consider  the  case  when  the  drain  is  buried  in  a  stratified  soil  at  an  interface 
as  shown  in  figure  42,  a.  The  upper  half  of  the  drain  is  in  a  layer  of  conduc- 
tivity Kx  and  the  bottom  half  of  the  drain  is  in  a  layer  of  conductivity  Kg.  The 
problem"  can  be  solved  with  the  same  procedure  described  above.  The  only 
difference  is  the  addition  of  the  term  Cd  which  is  reported  in  table  1.  The  net- 
work representing  each  block  of  soil  is  shown  in  figure  42,  b,  where  m  is  the 
ratio  Kj/Kb.  The  voltage  to  be  determined  is  Vo,  which  is  shown  in  the  center 
of  figure  42,  c.  Using  the  same  procedure  described  above,  the  voltage  Vo  can 
be  determined  as   shown  in  equation  53. 

.   2mCdV|  »  Cd(m  +  l)V2  +  2CdV3    +  2(m  +  I  )V4  } 
-°  "                               (m+  l)(3Cd  +2) 


"Rectangular  blocks  of  soil  around  the  drain  region  can  be  handled  by  proper  manipulation  of  the  equation  in 
(12,  p.  593,  type  of  line  T).    This  is  done  by  dividing  Z0  by  2  and  changing  the  term  JJ  toW/2. 
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In  the    case  when  the    drain   is  buried  within  a  layer    so  that  m  =  1,  the  expres 
sion  for  the  potential  adjacent  to  the  drain  is  given  by  equation  54. 


Yo  " 


v,  +  Y2  +  y3 


+  £fr 


»♦* 


(54) 


This  equation  is  of  similar  form  to  Luthin's  equation  4  in    reference  (8),  where 

the  term   h   in  Luthin's  equation  4    is  equivalent  to    our  ^S-  .    However,  Luthin's 

expression  is  based  on  a  linear  change  in  the  voltage  between  the  drain  and  an 
adjacent  node;  whereas,  equations  53  and  54  assume  a  logarithmic  relation- 
ship. It  is  shown  in  another  report  (11)  that  a  logarithmic  relationship  is  nec- 
essary to  give  satisfactory  agreement  between  results  obtained  with  the  network 
and  analytical  solutions. 
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FIGURE  42.  — Diagram  showing  building  blocks  of  different  hydraulic  conductivity 
around  the  drain.  The  diagram  is  used  to  determine  the  potential  V  0  at  the 
nodes  nearest  the  drain. 
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m.     SELECTING  THE  CHARACTERISTIC  RESISTANCE  OF  THE  NETWORK 

There  are  two  factors  which  must  be  considered  in  choosing  the  charac- 
teristic resistance  R0  of  the  network:  These  factors  are  concerned  with  (1) 
the  error  introduced  in  voltage  measurements  and  (2)  the  error  resulting  from 
contact  resistance  in  the  many  jack  and  plug  assemblies. 

In  voltage  measurements  a  voltmeter  is  connected  in  parallel  with  the 
circuit.  Since  a  voltmeter  has  an  internal  resistance  of  approximately  11 
megohms,  connecting  the  voltmeter  to  a  circuit  having  a  resistance  RQ  will 
yield  a  total  resistance  R  as  follows  (see  equation  6): 

..  ll,000,000Ro 
-  "  11,000,000+Ro 

This  will  introduce    an   error  in  the  measurement  of   R0  of    -r   —        100  percent, 

100  R  ~  — ° 

or    .i    MnTT^imm      percent.     The    latter    expression  indicates  that  if  R0  is  large, 
n,uuu,uuU+R0  —  6 

one   will  have    a   high  percentage    of   error.     It    is   found  that  when  a   network  of 

resistors    is    assembled   to    represent   a    soil  profile,    the    overall  resistance  of 

the  network  is  in  the    neighborhood  of  RQ.     It  follows,  therefore,  that  R0  should 

be  chosen  small  to  reduce  errors  in  voltage  measurement. 

The  second  factors  concerns  errors  resulting  from  wiring  resistances  and 
collective  resistances  in  plug  and  jack  assemblies.  The  overall  resistance 
added  by  these  sources  is  about  1  to  5  ohms.  If  the  network  resistance  is  around 
Ro  ohms  and  the  average  wiring  resistance  is,  say,  2.5  ohms,  the  percentage 
error  due  to  wiring  resistance  is  as  follows: 


/(Ro  +  2.5)  -  Rp\x  1 

I     so     ; 


00  =  250  percent 
Ro 


From   this    expression,    one    can  see    that    if  Ro  is    small,   the  percentage  error 
will  also  be  large. 

With  the  two  sources  of  error  in  mind,  we  can  determine  an  optimum  value 
for  the  characteristic  resistance  R0  of  the  network.  This  is  done  by  first 
averaging  the  values  of  the  logarithm  of  1 1  megohms  (internal  resistance  of 
the  voltmeter)  and  of  2.5  ohms  (wiring  resistance).  Then  one  takes  the  antilog 
of  the  average  to  obtain  the  value  of  the  characteristic  resistance  RQ.     That  is, 

log  11,000,000  =  7  and  log  2.5    «  0.4  so   that  the  antilog  of  7  +2°A  =  5,000  ohms. 

If  an  R0  of  5,000  ohms  is  chosen,  the  percentage  error  due  to  internal 
resistance  of  the  voltmeter  (upper  limit)  is  0.045  percent  and  the  percentage 
error  due  to  wiring  resistance  (lower  limit)  is  0.050  percent.  Any  deviation 
of  R0  from  5,000  ohms  will  increase  one  and  decrease  the  other.  For  example, 
if  an  R0  of  200  ohms  is  chosen,  the  error  of  the  upper  limit  is  0.00018  percent 
and  the  error  of  the  lower  limit  is  increased  to   1  percent. 
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Using  a  Voltmeter  to  Measure  Currents 


One  can  use  a  digital  voltmeter  to  measure  current  flow.  This  is  done  by- 
attaching  a  standard  resistor  of  1,10,  or  100  ohms  in  series  with  the  circuit 
having  resistance  of  Ro  ohms  and  then  using  the  voltmeter  to  measure  the  volt- 
age across  the  standard  resistor  as  shown  in  figure  43.  The  magnitude  of  R  is 
not  important,  just  as  long  as  the  voltage  across  it  is  the  voltage  desired.  The 
current  going  through  both  the  standard  resistor  and  the  resistance  R  will  be 
the  same  since  they  are  in  series.  By  reading  the  voltage  drop  across  a  stand- 
ard resistor,  the  current  I  can  be  determined  as  — ==. 


R  standard 
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FIGURE  43.  — Diagram  showing  the  use  of  VTVM  to  measure  current  flow. 


Error  in  Setting  Variable  Resistors  with  a  Wheatstone  Bridge 


In  using  the  Wheatstone  bridge  to  set  a  resistor  to  the  desired  value,  one 
proceeds  as  follows:  In  two  arms  of  the  bridge  place  precision  resistors  (0.05 
percent)  of  equal  magnitude.  Place  resistance  decade  box  of  0.05  percent 
accuracy  covering  a  range  from  0.1  ohm  to  1  megohm  in  another  arm  of  the 
bridge.  Then,  insert  the  resistor  to  be  set  into  the  remaining  arm  (fig.  44).  To 
set  the  resistance,  set  the  decade  box  to  the  desired  value,  and  adjust  the  un- 
known resistor  until  the  voltage  V  is  zero.  However.it  is  impractical  to  adjust 
the  unknown  resistor  until  the  voltage  A  V  reads  exactly  zero.  In  such  case,  the 
desired  resistance  will  have  an  error  of  A  R.  To  find  the  percentage  error,  one 
proceeds  as  follows:     Let  the  two  fixed  resistances  and  the  decade  box  all  have 
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FIGURE  44.— Diagram  showing  the  current  of  the  Wheatstone  bridge.  The 
unknown  resistor  which  is  to  be  set  to  a  certain  value  is  put  in  the  arm 
marked  "unknown"  and  is  adjusted  until  the  voltage  AV  approaches  zero. 


the  same  resistance  R.    Apply  a  voltage  V  to  the  circuit.     The  amount  of  error 
he  resistance  can  b 
AR  x  100      400AV 


aR      4  AV 

in  setting  the  resistance  can  be  expressed  as^=-  =  — v~>  so  that  the  percentage 


error  is 


R 


V 


In  setting    each    resistance,    apply    a   voltage  V  of  not   more   than  30  volts  to  the 

bridge    and  keep   the    residual  AV   below    0.025   volt.     The   percentage    error    in 

.  .    .  .  .  .        .      .         400  x  0.025        A  ^^ 

setting  each  resistance  can  be  estimated  as  ^ =  0.33  percent. 

In  the  network    system,    these    resistors  are   joined   in   mesh  and  the    maximum 

S  A  R 

errors    should   not   exceed   the    amount =  where    n   is   the    number    of   the  re- 

H  — 

sistors  and  A  g  is  the  error  of  each  resistor. 
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